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ARTICLE INFO ABSTRACT 


Keywords: Although canine respiratory coronavirus (CRCoV) is an important respiratory pathogen that is prevalent in many 
Betacoronavirus countries, only one complete genome sequence of CRCoV (South Korea strain K37) has been obtained to date. 
Canine respiratory coronavirus Genome-wide analyses and recombination have rarely been conducted, as small numbers of samples and limited 
ooneme genomic characterization have previously prevented further analyses. Herein, we report a unique CRCoV strain, 
Origin . . denoted strain BJ232, derived from a CRCoV-positive dog with a mild respiratory infection. Phylogenetic 
Phylogenetic analysis é £ : : 7 

Recombination analysis based on complete genome of all available coronaviruses consistently show that CRCoV BJ232 is most 


closely related to human coronavirus OC43 (HCoV-OC43) and BCoV, forming a separate clade that split off early 
from other Betacoronavirus 1. Based on the phylogenetic and SimPlot analysis we propose that CRCoV-K37 was 
derived from genetic recombination between CRCoV-BJ232 and BCoV. In detail, spike (S) gene of CRCoV-K37 
clustered with CRCoV-BJ232. However orflab, membrane (M) and nucleocapsid (N) genes were more related to 
Bovine coronavirus (BCoV) than CRCoV-B232. Molecular epidemic analysis confirmed the prevalence of CRCoV- 
BJ232 lineage around the world for a long time. Recombinant events among Betacoronavirus 1 may have 
implications for CRCoV transmissibility. All these findings provide further information regarding the origin of 
CRCoV. 


1. Introduction CRCoV was included in species Betacoronavirus 1, which also 


contain BCoV, equine coronavirus, HCoV-OC43, porcine hemaggluti- 


Canine respiratory coronavirus (CRCoV), an enveloped positive- 
stranded RNA virus, belongs to the betacoronavirus genus, subgroup 2a 
(Erles et al., 2007; Erles et al., 2003). Unlike the enteric canine 
coronavirus (CCoV) belonging to the alphacoronavirus genus, CRCoV 
is associated with mild to severe respiratory signs and proposed as an 
etiological agent of canine infectious respiratory disease (CIRD) 
(Buonavoglia and Martella, 2007). CRCoV was first detected in 2003 
by RT-PCR in the UK from dogs suffering from severe respiratory 
disease (Erles et al., 2003). Later, the CRCoV-4182 strain was success- 
fully isolated from a dog with respiratory signs and was cultured in the 
HRT-18 cell line (Erles et al., 2007). To date, data have been made 
publicly available for only one complete genomic sequence (CRCoV- 
K37 strain), which was isolated from a South Korea dog (An et al., 
2010; Lim et al., 2013). 


nating encephalomyelitisvirus (PHEV), human enteric coronavirus 
(HECoV) and CRCoV (King, 2011). In contract, Canine coronavirus 
(CCoV) is a member of alphacoronavirus, which include feline corona- 
virus and porcine transmissible gastroenteritis among others. The virus 
was distinct from CCoV and only showed 69% nucleotide identity in 
polymerase and 21% amino acid identity in spike protein (Buonavoglia 
and Martella, 2007; Erles and Brownlie, 2008). 

The evolutionary origin of CRCoV remains uncertain, although 
successful experimental infection of puppies with BCoV has suggested 
that CRCoV is of bovine origin and transmitted to dogs from cattle 
(Kaneshima et al., 2007). Meanwhile, HCoV-OC43 may have emerged 
after viral transmission from cattle to people (Vijgen et al., 2005). The 
high nucleotide identity between CRCoV and BCoV, HCoV-OC43 and 
HECoV indicated that the four viruses have a common ancestor (Erles 
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and Brownlie, 2008), and demonstrated the occurrence of repeated 
host-species shifts (Vijgen et al., 2006). Although CRCoV strains are 
frequently circulating in farmed and pet dogs worldwide (Erles and 
Brownlie, 2005; Kaneshima et al., 2006; Priestnall et al., 2006; 
Priestnall et al., 2007; Yachi and Mochizuki, 2006), causing respiratory 
disease, little is known about their genomic evolution and recombina- 
tion. In the present study, we sequenced and analyzed the first complete 
genome of CRCoV derived from China (CRCoV-BJ232 strain) and 
discuss its role in CRCoV evolution. 


2. Materials and methods 
2.1. Ethics statement 


This research was reviewed and approved by the ethics committee 
of the National Institute for Viral Disease Control and Prevention of the 
Chinese CDC. All animals were treated strictly according to the guide- 
lines for the Laboratory Animal Use and Care from the Chinese CDC. 


2.2. Specimen collection and handling 


From December 2013-March 2014, 246 swab specimens (nasophar- 
yngeal and anal swabs) were collected from diseased dogs at the Animal 
Hospital of China Agricultural University. Most of the dogs showed 
symptoms of gastrointestinal and/or respiratory tract infections, in- 
cluding coughing, lacrimation, vomiting and diarrhea. All the swab 
samples were collected by professional veterinarians and were stored at 
— 80 °C. 


2.3. CRCoV detection 


For molecular detection, total RNA was extracted from swab speci- 
mens using a QIAamp Viral RNA Mini Kit (QIAGEN, Germany) 
according to the manufacturer’s instructions. The RNA was used as a 
template for reverse transcription-PCR (RT-PCR), and the partial 
sequences of the S, RNA-dependent RNA polymerase (RdRp), and N 
genes were amplified using gene-specific primers (Table $1). In detail, 
reverse transcription was performed by using a SuperScript III kit 
(Invitrogen, San Diego, CA). The PCR mixture (25 ul) contain 1 pl 
cDNA, 1 ul forward primer(10 UM), 1 pl reverse primer (10 LM), 12.5 pl 
Q5° High-Fidelity 2X Master Mix (New England Biolabs) and 9.5 ul 
nuclease-free water; the mixtures were amplified with 30 cycles at 98 °C 
for 10s, 52 °C for 30s, and 72 °C for 1 min and final extension step at 
72°C for 2min in an automated thermal cycler (BIO-RAD, T100 
Thermal Cycler). 


2.4. Complete genome sequencing of CRCoV 


To better characterize the CRCoV discovered here, the CRCoV- 
positive samples were pretreated and used for random amplification as 
described previously (Wang et al., 2016); the amplified DNA was used 
as a template for Illumina HiSeq2500 sequencing (2 x 125 bp, paired- 
end reads). The sequencing procedure was conducted by Beijing Berry 
Genomic Company. The full-length genome of CRCoV was assembled 
using SOAP de novo and was annotated based on the sole annotation of 
CRCoV K37 strain (JX860640.1). 


2.5. Phylogenetic analysis 


To determine the evolutionary relationships between the CRCoV 
discovered here and those identified previously, we collected most of 
the complete sequences of the Betacoronavirus 1 and conducted 
nucleotide acid sequence alignments using the MAFFT algorithm. The 
phylogenetic analyses were conducted based on the nucleotide acid 
sequences of the orflab, S, M, and N genes and complete genome of 
CRCoV. All the phylogenetic trees were constructed in the PhyML 
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software using the maximum likelihood method, with 500 bootstrap 
replicates. The following Betacoronavirus 1 genome sequences from the 
GenBank database were used for the phylogenetic analysis as follows: 
BCoV reference strain (NC_003045), HECoV (NC_012950), HCoV-OC43 
(NC_005147), CRCoV-37 (JX860640), and BCoV strain Quebec 
(AF220295). 


2.6. Recombinant analysis 


To confirm the potential recombination events, similarity plots and 
bootscanning analyses were conducted using the SimPlot software; a 
sliding window of 1000 nucleotides and 100-nucleotide steps were used 
as the default settings. In addition, a schematic diagram of the complete 
CRCoV genome was examined to determine the recombination sites. 


2.7. Nucleotide sequence accession number 


The complete genome of ChinaBJ232 was deposited in the GenBank 
database under accession no. KX432213. Other sequences generated in 
this study, including the partial S, RdRp and complete N gene 
sequences, were also assigned accession numbers: KX432214- 
KX432217. 


3. Results 
3.1. Identification of CRCoV 


A total of 246 swab samples were collected in Beijing from 
December 2013 to March 2014.RT-PCR was performed to detect 
CRCoV, and 16 (6.5%) out of 246 swabs originating from Chinese dogs 
tested positive for CRCoV. In addition, five CRCoV-positive samples 
with high titers were initially processed for complete genome sequen- 
cing using deep sequencing technology. 


3.2. Complete genomic characterization of the BJ232 strain 


Due to the limited specimen size and low viral load, five clinical 
samples were used for deep sequencing, and only one complete 
genomic sequence was eventually obtained. The genomic data, formally 
referred to as “Canine Respiratory Coronavirus strain BJ232” and in 
this report as “CRCoV-BJ232,” were deposited in GenBank (KX432213). 
The CRCoV-BJ232 genome isolated in China was found to have 30,868 
nucleotides and showed a genomic organization that was highly similar 
to the Betacoronavirus l1species: 5’-UTR(nt 1-143),ORFlab (nt 
144-21427), ORF2 (nt 21437-22273), HE gene (nt 22285-23559), S 
gene (nt 23574-27665),E gene (nt 28355-28609), M gene (nt 
28624-29316), N gene (nt 29326-30672), and 3’-UTR (nt 
30673-30868). Between the S and E gene, three small ORFs were 
identified: ORF4a (nt 27655-27789), ORF4b (nt 27825-27902) and 
ORF4c (nt 28039-28293), encoding 4.9-, 2.7- and 12.8-kDa non- 
structural proteins, respectively, which were analogous to those noted 
in strain CRCoV-K37. Meanwhile, conventional PCR and Sanger 
sequencing were conducted to verify the complete genome sequence 
of CRCoV-BJ232 with 100% nucleotide identity. Thus far, there are 
only two CRCoV strains with full-length genomic information submitted 
to GenBank; the previously reported CRCoV-K37 and our current strain, 
CRCoV-BJ232. Unfortunately, attempts to isolate CRCoV-BJ232 from 
the five positive samples were unsuccessful; future studies to isolate the 
virus from more canine samples will allow characterization of its 
pathogenicity and transmissibility. 


3.3. Phylogenetic analysis 
Based on the analysis of whole-genome phylogeny, CRCoV-BJ232 


was classified into Betacoronavirus 1 species of betacoronavirus group 
2a (Fig. 1A) and shares a common origin with BCoV, CRCoV-K37, 
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NC 022643.1 Betacoronavirus Erinaceus/VMC/DEU/2012 isolate ErinaceusCoV/2012-216/GER/2012 


100-— NC 009019.1 Bat coronavirus HKUS-1 
100 f_ NC 008315. 1 Bat coronavirus (BtCoV/133/2005) 
NC028918.1 BtTp-BetaCoV/Gx2012 
100 NC 019843.3 Middle East respiratory syndrome coronavirus 
NC 028909.1 BtPa-BetaCov/GD2013, 
100 NC 009020.1 Bat coronavirus HKUS-1 
NC 006577.2 Human coronavirus HKU1 
NC 026011. 1 Betacoronavirus HKU24 strain HKU24-RO5005I 
100 100 p— AC 000192.1 Murine hepatitis virus strain JHM 
ae: NC 001846.1 Murine hepatitis virus strain A59 
NC 012936.1 Rat coronavirus Parker 
100 NC.017083,1 Rabbit coronavirus HKU14. 
100) NC 010327.1 Equine coronavirus 2a 
NC 007732.1 Porcine hemagglutinating encephalomyelitis virus 
NC 005147.1 Human coronavirus OC43 
© Canine respiratory coronavirus strain 8J232 
NC 003045. 1 Bovine coronavirus 
NC 012950.1 Human enteric coronavirus strain 4408 
1+ JX860640 Canine respiratory coronavirus strain K37 
100 NC 001451.1 Avian infectious bronchitis virus 
NC.010800.1 Turkey coronavirus 
NC 010646.1 Beluga Whale coronavirus SW1 
100 NC 016994.1 Night-heron coronavirus HKU19 
100 NC 016995.1 Wigeon coronavirus HKU20 
NC 016996.1 Common-moorhen coronavirus HKU21 
a NC 011549. 1 Thrush coronavirus HKU12-600 
NC 0169911 White-eye coronavirus HKU16 5 
La NC.016993.1 Magpie-robin coronavirus HKU18 
NC011550.1 Munia coronavirus HKU13-3514 
NC 016992. 1 Sparrow coronavirus HKU17 
100 NC 016990.1 Porcine coronavirus HKU15 
100 NC 028824.1 BtRf-AlphaCoV/YN2012 
NC 009988.1 Bat coronavirus HKU2 
100 NC 028806.1 Swine enteric coronavirus strain Italy/213306/2009 
NC 002306. Feline infectious peritonitis virus 
100 NC 023760.1 Mink coronavirus strain WD1127 
100 NC 028752 1 Camel alphacoronavirus isolate camel/Riyadh/Ry141/2015 
00, NC 002645.1 Human coronavirus 229E 
59 NC 005831.2 Human coronavirus NL63 
100 NC 009657.1 Scotophilus bat coronavirus 512 
NC 003436.1 Porcine epidemic diarrhea virus 
100 NC 022103.1 Bat coronavirus CDPHE 15/USA/2006 
NC 028811.1 BtMr-AlphaCoV/SAX2011 
100 p— NC 028839.1 BtMs-AlphaCoV/GS2013 a 
100 NC 028814.1 BtRf-AlphaCoV/HuB2013 
NC 018871.1 Rousettus bat coronavirus HKU10 
NC.028833.1 BtNv-AlphaCoV/SC2013 
98 98-— NC 028756.1 BIMf-AlphaCoV/AH2011 
100 |_ NC 010437.1 Bat coronavirus 1A 
100 NC028761.1 BMf-AlphaCoV/JX2012 
NC 010436.1 Bat coronavirus 1B 
100 p>— NC 028771.1 BtMt-AlphaCoV/GD2012 
NC 010438.1 Bat coronavirus HKUB 
oa NC 028781.1 BIMf-AlphaCoV/HuB2013 
NC 028797. 1 BIMf-AlphaCoV/HeN2013 
100 + NC 028790.1 BtMI-AlphaCoV/FJ2012 
NC 022787.1 Porcine torovirus strain SH1 
100 NC 007447.1 Breda virus 
NC 026812.1 Chinook salmon bafinivirus isolate NIDO 
100 NC 008516.1 White bream virus 


at 
O14 


2c 


NC 006577.2 Human coronavirus HKU1 


NC 026011.1 Betacoronavirus HKU24 strain HKU24-R050051 
AC 000192.1 Murine hepatitis virus strain JHM 

NC 001846.1 Murine hepatitis virus strain A59 

NC 012936.1 Rat coronavirus Parker 

NC 017083.1 Rabbit coronavirus HKU14 

NC 010327.1 Equine coronavirus 

NC 007732.1 Porcine hemagglutinating encephalomyelitis virus 

NC 005147.1 Human coronavirus OC43, 

@ X432213 Canine respiratory coronavirus strain BJ232 Pane 
NC.003045.1 Bovine coronavirus 

97 NC 012950.1 Human enteric coronavirus strain 4408 


Sil JX860640 Canine respiratory coronavirus strain K37 


0.05 


Fig. 1. Phylogenetic analysis of all available coronaviruses with reference sequence based on complete genome. A, the current taxonomy of the family Coronaviridae; B, the evolution of 
Betacoronavirus 1 based on complete genomes. StrainBJ232 (accession no. KX432213, indicated by the red dot) was the first CRCoV with complete genome derived from China. The 
previously reported CRCoV-K37 strain, which was detected in South Korea, was the sole other CRCoV strain with complete genome sequence available. Some other representative strains 


of the Betacoronavirus genus, subgroup 2a, were also analyzed in PhyML program using the maximum likelihood method, with 500 bootstrap replicates. The numbers at the nodes 


represent the bootstrap support. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 


HECoV and HCoV-OC43. Moreover, the latter four viruses clustered 
together with high branch bootstrap values (Fig. 1B). To investigate the 
genetic relationships among the Betacoronavirus 1 species, we gener- 
ated phylogenetic trees based on the sequences oforflab, S, M and N 
genes (Fig. 2). Phylogenetic analysis based on different genes consis- 
tently show that CRCoV BJ232 is most closely related to HCoV-0C43, 
forming a separate clade of CRCoV that split off early from other 
Betacoronavirus 1. However, CRCoV-K37 strain closely clustered with 
CRCoV-BJ232 strain in S gene, and was more related to BCoV strain 
Quebec (AF220295.1) in the region of orflab, M and N genes. The 
major genes of the CRCoV-BJ232 strain showed a consistent trend in 
the evolution, while the phylogenetic topology conflicted in the sub- 
regions of CRCoV strain K37, especially in the S and orflab genes. 
These findings are consistent with recombination, a phenomenon not 
uncommon in coronaviruses, implying that the CRCoV K37 strain may 
have arisen from genetic recombination. 


3.4. Recombinant analysis 


Similarity plots and bootscanning analyses (Fig. 3) were performed 
to confirm the recombination between CRCoV and BCoV. The pre- 
viously reported CRCoV K37 strain was used as a query sequence, and 
was compared with HCoV-OC43 (NC.005147), BCoV (BCoV- 
RefNC_003045, BCoV-DB2 DQ811784, BCoV-Quebec AF220295), HE- 
CoV (HECoV-Ref NC_012950, HECoV-4408 FJ415324), and CRCoV 
(CRCoV-BJ232 KX432213). In the S gene, strain K37 possessed the 
highest similarity with the BJ232 strain between position 23,842 and 
29,072, as expected. However, in the orflab, HE, M and N genes, strain 


K37 was apparently related to BCoV, which further confirmed the 
occurrence of genetic recombination between CRCoV and BCoV. The 
consistency in the results of bootscanning and phylogenetic analyses 
also supported the possibility of recombination in the S gene. 


3.5. Prevalence of recombinant CRCoV worldwide 


We collected all available CRCoV sequences from the GenBank 
database and conducted comprehensive phylogenetic analyses sepa- 
rately based on HE, S, M and N genes to understand the epidemiology 
distribution of the recombinant strain(South Korea strain K37). The 
BJ232, BJ221, BJ34 and BJ202 strains were collected from China in 
this study, while the K37, K39 and K9 strains came from South Korea, 
and the 240/05 and 4182 strains were derived from Italy and the 
United Kingdom, respectively. As shown in Fig. 4, the CRCoV strains 
from China, including the BJ232, BJ221, BJ34 and BJ202 strains, 
clustered with strains prevalent in Italy and the United Kingdom (the 
240/05 and 4182 strains), and the CRCoV strains detected in South 
Korea (K37, K39 and K9 strains) were all derived from recombinant 
lineage. 


4. Discussion 


CRCovV is a causative agent for CIRD and a novel pathogen detected 
in respiratory samples from dogs (Erles and Brownlie, 2008). A 
serological survey of CRCoV conducted worldwide revealed that 
antibodies against CRCoV were common in the canine population 
(Erles et al., 2004; Isaacs et al., 1983; Pratelli et al., 2003; Priestnall 
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Fig. 2. Evolutionary analyses of major coding genes of Betacoronavirus 1. The coding genes used in the phylogenetic analyses included orflab, S, M and N. The two CRCOV strains (BJ232 
and K37 strains) are indicated in red. All the sequences were analyzed in PhyML program using the maximum likelihood method, with 500 bootstrap replicates. 


et al., 2006). In contrast, corresponding genome-wide analyses have 
rarely been conducted (An et al., 2010; Lorusso et al., 2009). In the 
present study, we characterized the first complete genome of CRCoV 
(named strain BJ232) derived from China, and phylogenetic analyses 
showed that CRCoV-BJ232occupied a deep branch at the root of 
members of Betacoronavirus 1, being distinct from HCoV-OC43 and 
BCoV. Another previously described CRCoV strain, K37, was identified 
to have arisen from recombination between CRCoV-BJ232 and BCoV. In 
detail, the CRCoV-K37 strain was related to CRCoV-BJ232 in the S gene, 
while in the orflab, HE, M, and N genes, CRCoV-K37 was noted to 
possess a close relationship with BCoV. 

Recombination is a common phenomenon in coronaviruses and 
thought to contribute to the emergence of new pathotypes (Gorbalenya, 
2008; Wang et al., 2015). Thus far, most of the recombination events of 
coronaviruses have been reported between species of the same group 
(Herrewegh et al., 1998; Keck et al., 1988), such as among the bat- 
associated CoVs (Corman et al., 2014) and 229E-related CoVs (Corman 
et al., 2015), and the major recombination breakpoint has mainly been 
within the S gene. Genetic recombinant among alphacoronavirus 1 
species frequently occurred, such as ferret coronaviruses (Lamers et al., 


2016; Minami et al., 2016). However, few researches about genetic 
recombination of Betacoronavirus 1 species have been conducted. The 
discovery of genetic recombination between CRCoV and BCoV confirms 
the phenomenon, and intraspecies recombination therefore exists for 
Betacoronavirus 1 species. 

Phylogenetic analyses based on all available CRCoV genes showed 
that the two clusters of CRCoV, including recombinant strain (South 
Korea K37 strain) and the Chinese BJ232 strain, have been prevalent 
worldwide for many years. CRCoV and BCoV were closely related in 
Betacoronavirus 1, located on the outer leaves of phylogenetic tree, and 
form an independent branch. It was difficult to determine the evolution 
orientation (transmission from cattle to canine or both from murine), 
although the latest research supports the murine origins of 
Betacoronavirus 1 (Lau et al., 2011, 2015). More CRCoV genome 
sequences are needed to understand the evolutionary relationship 
between CRCoV and BCoV. 

Regarding the pathogenicity and transmissibility of CRCoV, the 
Chinese BJ232 strain and South Korea K37 strain were both derived 
from diseased dogs with respiratory infections. Genomic analyses 
showed that the S genes, which mediate coronavirus attachment, were 
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Fig. 3. Genetic recombination analyses of the complete genome of the South Korea K37 strain. (A) A schematic diagram of the complete genome of the CRCoV-K37 strain. (B) The results 
of a SimPlot similarity analysis. (C) The results of a bootscanning analysis. The CRCoV-K37 strain was used as the query sequence and was compared with seven other representative 
strains of BCoV-like viruses, including HCoV-OC43 (NC_005147.1), BCoV (NC_003045.1), BCoV DB2 strain (DQ811784), HECoV (NC_012950.1), HECoV 4408 strain (FJ415324.1), BCoV 
Quebec strain (AF220295.1), and CRCoV-BJ232 (KX432213). A window size of 1000, a step size of 1000 bp and 100 replicates were used as the default settings. 
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Fig. 4. Phylogenetic analyses of all available CRCoV strains based on HE, M, N and S genes. The BJ232, BJ221, BJ34 and BJ202 strains come from China in this research,K37, K39 and K9 
strains were collected from South Korea, and the 240/05 and 4182 strains were derived from Italy and the United Kingdom, respectively. All the genes were analyzed by PhyML program 


using the maximum likelihood method, with 500 bootstrap replicates. 


highly similar in these strains. Because recombination may be a driving 
force for the formation of pathogenic viruses from less pathogenic virus 
(Forni et al., 2017), recombinant may play an important role in the 
pathogenicity and transmissibility of CRCoV. Limitation of this study 
was that the isolation of CRCoV-BJ232 failed on cell culture. More 
CRCoV strains were needed to be isolated or rescued using reverse 
genetic technique to perform experimental inoculation of dogs. 

In conclusion, herein, we reported the first complete genome 
sequence of a CRCoV derived from a diseased dog in China, and 
confirm that intraspecies recombinant exists for Betacoronavirus 1 
species. The results may have implications for the origin and evolution 
of CRCoV. 
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